Recent results on anomalous J/iJj suppression in Pb-Pb collisions at 158 AGeV/c 
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We report our latest results on charmonium suppression as measured at the CERN-SPS in 
Pb-Pb interactions at an incident beam momentum of 158 GeV/c per nucleon. Preliminary 
results obtained from the most recent sample of data collected in year 2000 under improved 
experimental conditions are compared with published results. For the most peripheral Pb- 
Pb collisions, J/tp suppression agrees with the normal absorption measured from interactions 
of lighter nuclei while a steady increasing abnormal supression is observed with increasing 
centrality. 



1 Introduction 



The study of charmonium production in ultrarelativistic heavy ion collisions by experiments 
NA38 and NA50 has been motivated by the search for a phase transition of nuclear matter from 
its normal state to a deconfined quark gluon plasma, as predicted to occur, under extreme energy 
densities or temperature conditions, by non-perturbative QCD. In such a deconfined matter, the 
J/ijj is expected to be suppressed by Debye colour screening -'^ . As the vector- meson is produced 
at the very early stages of the collision, its leptonic decay provides a probe of these stages, since 
this mode is insensitive to strong interactions and, therefore, to the hadronization stage. 

2 Data selection and analysis method 

The main element of the NA50 detector is a muon spectrometer. It is complemented by central- 
ity detectors which measure the neutral transverse energy Et-, the beam spectators energy Ezdc 
and the charged multiplicity produced in the interaction. Several beam monitoring devices are 
also used for various purposes. A detailed description can be found inl^. The goal of the year 
2000 Pb-Pb data collection is to investigate, in particular, the most peripheral collisions which, 
in previous data samples, could be contaminated by Pb-air interactions. The data have been 
taken, therefore, with the target placed in the vacuum pipe. 

The event selection is done according to the following criteria: dimuons are selected in the rapid- 
ity window 2.92 < yiab < 3.92 (0 < ycM < 1) and with a Collins-Soper angle | cosOcs I < 0.5, 
leading to an acceptance of around 14% in the mass region of interest. Parasitic interactions of 
an incident Pb-ion in the Beam Hodoscope and upstream from the target are rejected by a BH 
interaction detector and by Anti-halo counters. On-target interactions are selected requiring a 
good correlation between hits in the two planes of the Multiplicity Detector. Double interactions 
are discarded by a shape analysis of the signal in the Electromagnetic Calorimeter. Residual 
pileup events are rejected by a 2a cut in the Et — Ezdc correlation!^. Off-target muons are 
excluded by a cut on the transverse distance between the extrapolated muon track and the beam 
line. 

In each centrality bin of the analysis, the opposite-sign dimuon mass spectrum is adjusted in 
the pertinent mass region, accounting for the different ingredients: J/il^, ip' and the continuum 
formed by Drell-Yan, DD and the combinatorial background from vr and K decays. The shapes 
of the resonances are obtained by Monte Carlo simulation followed by reconstruction with the 
same program as used for the real data analysis. Drell-Yan and DD are generated at leading 
order by Pythia and treated by the NA50 simulation-reconstruction chain. The combinatorial 
background is properly extracted from like-sign pairs through the formula Nbg = 2\/ N^^N . 

3 Before the year 2000 The state of the art 

The NA38/NA51/NA50 experiments have verified that, as expected theoretically, the Drell-Yan 
cross-section cr(DY) is proportional to the number of elementary nucleon-nucleon collisions, from 
proton-proton up to lead-lead interactions included (see Fig. The ratio of the cross-sections 
for J/'ijj and Drell-Yan production is therefore proportional to the J/ip production cross-section 
per nucleon-nucleon collision. Moreover, the ratio B^^a{J / (t{DY) is directly measured in 
the so-called standard analysis method and insensitive to reconstruction or trigger inefficiencies 
(identical muon pair topology for J/ij: and Drell-Yan events) or to absolute normalization un- 
certainties which, obviously, cancel out. Fig.|21shows the ratio B ^aoiJ / 'ip) / oiDY) as a function 
of Et, for the previously collected Pb-Pb data samples - 1995 M 1995 El and 1998^1. The 
closed points are the results obtained with the standard analysis method described above while 
the open points are the 1996 and 1998 results obtained through the "Minimum Bias analy- 
sis" method explained in [4,5]. The continuous line is the J/ if) normal absorption in nuclear 



matter and fits the NA51 p-p and p-d and NA38 p-A and S-U results. A clear departure of 
the B ^^a{J / ip) / a{DY) ratios from this normal absorption curve is observed at mid-centrality, 
suggesting the onset of an anomalous J /ip suppression mechanism. Futhermore, no saturation 
of the ratio at high Et is observed, as opposed to the absorption curve trend. 
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Figure 1: Drell-Yan K-factor as a function 
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4 The preliminary year 2000 results 



» PbPb 1995 
■ PbPb 1996 
• PbPb 1998 

□ PbPb 1996 with Minimum Bias 
o PbPb 1998 with Minimum Bias- 
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Figure 2; J/tp/DY as a function of Et; the 
absorption curve is a fit to the NA38 p-A 
and S-U data 



The NA50 Collaboration has a new high statistics p-A data sample, obtained with a 450 GeV/c 
incident energy beam. Since the J/^l^ absorption cross sections for p-A and S-U collisions are 
compatible within the statistical errors, a simultaneous Glauber fit is performed which leads 
to the common value of (j{J /ip)abs = 4.3 it 0.5 mb This result, combined with the NA38 
S-U data reanalyzed in order to have the same analysis conditions, leads to an updated normal 
nuclear absorption. 

From the above described fit to the dimuon invariant mass spectra, we extract the J/ij) and 
the Drell-Yan yields (the latter in the mass interval 2.9 - 4.5 GeV/c^) and obtain the ratio 
Bij,j_i(T{J /iIj) / a{DY) and its dependence as a function oi Et or Ezdc- Fig.Olshows the agreement 
among independent analyses with slightly different data selections and fit methods. It also shows 
that peripheral collisions indeed follow the absorption curve. On the other hand, the Pb-Pb 2000 
data confirm the departure from this normal behaviour of the J/ip/DY ratios at mid-centrality, 
already observed in the previously collected data. The results also confirm the steady decrease 
of the J/V' production rate at high Et- The same J/ip suppression pattern is seen in Fig. \^ 
when Ezdc is used as the centrality estimator. For this variable, the agreement is quite good 
between new and already published results, except for the most peripheral collisions, probably 
due to Pb-air interactions in the 1996 setup conditions. 

In the 2.9-4.5 GeV/c^ invariant mass range, muons from J /if) and Drell-Yan have similar 
features, namely acceptance and momentum. Nevertheless, different PDFs chosen to generate 
Drell-Yan lead to different mass distributions, as shown in Fig. El The difference under the J/ijj 
peak between the Drell-Yan cross sections generated with MRS43 and GRV LO amounts to 
10% systematic effect in the J j'^jDY normalizations, independent of centrality. To overcome 
this uncertainty, the 4.2-7.0 GeV/c^ mass range is chosen as the new J suppression reference. 
In fact, beyond 4.2 GeV/c^, the Drell-Yan pairs are the only contribution to the dimuon mass 
spectrum and thus both structure functions give the same normalization after adjustment to the 
data. 

Fig. El shows all NA50 results using the new mass range for Drell-Yan as reference. All data 
sets agree, regardless on the fact that Pb-Pb 1995 results were obtained using GRV LO PDFs 
while 1996, 1998 and 2000 used MRS43. Pb-Pb 1998 and 2000 analyses with GRV LO is also 
shown. 
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Figure 3: J/ip/DY as a function of Et 




O(abs)=4.3±0.5mb 



10 15 



20 



25 30 



35 



Figure 4: J/ip/DY as a function of Ezdc 




Figure 5: The Drell-Yan mass distribu- 
tions for MRS43 and GRV LO PDFs 
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Figure 6: J/ip/DY as a function of Et 



5 Conclusions 

In the most recently collected Pb-Pb data sample, the NA50 Collaboration has measured 
Bfj,fiCr{J/tp)/a{DY) as a function of centrality .With respect to previously collected data sam- 
ples, the use of a setup with the target region under vacuum allows to extend, confidently, the 
analysis to very peripheral collisions. 

These new data, confirming previous analyses results, show that J/ip production follows the 
normal nuclear absorption pattern for peripheral Pb-Pb collisions, while exhibiting a clear de- 
parture for central collisions and suggesting, therefore, the onset of another J ftp suppression 
mechanism. 
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